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In this work the connection between vortex condensation in a d-wave superconductor and 
'"^ , the QED3 gauge theory of the pseudogap is elucidated. The approach taken circumvents 

£^ ■ the use of the standard Franz-Tesanovic gauge transformation, borrowing ideas from 

Q ' the path-integral analysis of the Aharonov-Bohm problem. An essential feature of this 

O approach is that gauge-transformations which are prohibited on a particular multiply- 

connected manifold (e.g. a superconductor with vortices) can be successfully performed 
on the universal covering space associated with that manifold. 



1. Introduction 

Recently, much attention has been focused on understanding the pseudogap 1 - 2 < 3 ' 4 ' 5 



i phenomena of the high-temperature superconductors. One possible explana- 

Cn| [ tion 7 > 6 > 8 > 9 invokes the notion that the pseudogap is due to the presence of pairing 

correlations above the superconducting transition temperature T c . Within such a 
scenario, the lack of long-range phase coherence in the pseudogap is presumed to 
be due to the proliferation of vortex excitations in it. Following the approach of 
many recent papers S: 9 . 10 . 11 . 12 ^ we consider the problem of coupling vortices to the 
quasiparticles of a d-wave superconductor. As discussed in Refs. 10, 11, there is an 
q ' important distinction between the condensation of hc/2e (i.e. singly quantized) and 

hc/e (i.e. doubly quantized) vortices in a d-wave superconductor. Experimentally 13 , 
magnetic field-induced vortices seem to be exclusively of the singly quantized va- 
riety. We shall take this as evidence that, in considering vortex excitations in the 
pseudogap regime, it is sufficient to consider only hc/2e vortices. 

The technical issues associated with hc/2e vortices have been discussed in 
Refs. 10, 11, 14 and amount to the fact that certain singular gauge transforma- 
tions (which arise naturally when considering the Bogoliubov-de Gennes equation 
in the presence of vortices) lead to the presence of quasiparticle branch cuts. Here 
we shall adopt the point of view that these difficulties are not unique to the problem 
of performing gauge transformations in the presence of vortices: They arise when 
one considers the general problem of making gauge-transformations on multiply- 
connected manifolds. For example, consider the Aharonov-Bohm problem 15 , in 
which one imagines attaching current leads to a doubly-connected metallic ring 
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through which a solenoid penetrates. Even in the case of a solenoid with a radius 
sufficently small that the magnetic field is negligible in the ring, the presence of 
a nonzero vector potential leads to quantum-mechanical interference of electrons 
taking different paths around the ring while propagating between the leads. Thus, 
one cannot make a gauge-transformation which eliminates the effect of the vector 
potential on topologically distinct Feynman paths. 

In a similar sense, the winding of the pair-potential phase around a supercon- 
ducting vortex may not be trivially removed by a gauge transformation. In particu- 
lar, as noted above, making certain singular gauge-transformations in the presence 
of hc/2e vortices leads to quasiparticle branch cuts. Our aim here is to discuss a 
new approach to handling such branch cuts using the topological properties of path 
integrals. Recently, Franz and Tesanovic 14 have introduced a singular gauge trans- 
formation which avoids the introduction of quasiparticle branch cuts via a clever 
trick involving splitting the vortices into two distinct groups (A and B) and trans- 
forming the electrons relative to group A and the holes relative to group B. These 
authors find that the action governing the quasiparticle dynamics of a d-wave super- 
conductor in the presence of fluctuating vortices is given by the well-known problem 
of three-dimensional quantum electrodynamics (QED3): Dirac fcrmions coupled to 
a fluctuating "Berry" gauge field a. However, as we shall discuss, the propagator is 
not given by the usual propagator for QED3; it is given by the associated gauge- 
invariant propagator. This occurs because a is not a physical (i.e. electromagnetic) 
gauge-field; thus the propagator must not transform under gauge transformations 
on a. Gauge-invariant propagators have appeared many times in the context of 
strongly-correlated electron systems 16 > 17 > 18 ; hi particular the gauge-invariant prop- 
agator for QED3 exhibits intriguing non- Fermi liquid behavior n > 18 . 

To obtain further insight into the physics of fluctuating vortices in d-wave su- 
perconductors, here we shall take a novel approach inspired by Schulman's 19 topo- 
logical approach to the Aharonov-Bohm problem. Our task is to shed light on the 
connection between fluctuating vortices and fluctuating gauge-fields. Towards this 
end, we shall view a vortex as a "hole" which divides a superconductor into a 
doubly-connected space. As in the Aharonov-Bohm problem, the fact that the su- 
perconductor is multiply connected means that Feynman paths contributing to the 
quasiparticle propagator fall into topologically distinct sectors. By making different 
gauge transformations for topologically distinct Feynman paths, we shall see how 
the quasiparticle branch cuts may be expressed in terms of the Berry gauge field a. 
The purpose of this Paper is to arrive at the QED gauge theory of vortices in the 
pseudogap regime of the cuprates via a path integral technique. 

This Paper is organized as follows: In Sec. 2 we consider the problem of solv- 
ing the Bogoliubov-de Gcnnes (BdG) equation in the presence of a static array of 
vortices. Although our main interest will be in the problem of vortex fluctuations 
in a d-wave superconductor, the static case will be sufficient to motivate the tech- 
nical difficulties associated with vortices in superconductors. In Sec. 3, we review 
some results from the theory of the Aharonov-Bohm effect 15 from a point of view 
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due originally to Schulman 19 in which one considers path integrals in multiply- 
connected spaces. In Sec. 4, we revisit the BdG eigenproblcm from the point of 
view of gauge transformations in multiply-connected spaces and isolate the effect 
of the quasiparticle branch cuts on Feynman paths. In Sec. 5, these branch cuts 
are represented via a functional integral over an auxiliary field a; the theory is con- 
structed to be explicitly invariant under gauge-transformations associated with this 
field. In Sec. 6 we extend the results of Sec. 4 and Sec. 5 to the case of dynamic 
vortex excitations, finally arriving at an expression for the Bogoliubov-de Gcnnes 
propagator (in the presence of vortices) as a gauge-invariant Green function. In 
Sec. 7 we conclude with a brief discussion of our results. 



2. Bogoliubov-De Gennes Equation 

The QED 3 scenario of the pseudogap regime, like the nodal liquid scenario 10 which 
preceded it, focuses on the effect of vortex excitations on the quasiparticles of a 
d-wave superconductor. In the present section, we provide motivation by examining 
the problem of static vortex excitations in a superconductor. The quasiparticle 
excitations of a superconductor are described by the Bogoliubov-de Gennes (BdG) 
equation 

"(:rH(:;) <" 

/ - (V - teA) 2 - fx e W2Ae^/2 \ 

\ e-^Ae-W 2 (V + ieA) 2 + fi J ' { ' 

where fi is the chemical potential, ip is the local superconducting phase, e is the 
electronic charge, and A is the usual d-wave pairing operator, which we take to be 
given by A = A n p x p y . We have chosen units in which h 2 /2m = 1, with m being 
the quasiparticle mass. Here, u n and v n are the electron and hole parts of the BdG 
wavefunction, respectively. The electromagnetic gauge field A (which we shall often 
suppress) and ip conspire to give this theory the following local U(l) symmetry: 

Ap^Ap- (3) 

u n -> e lx u n , (4) 
v n -> e~ iX v n , (5) 
+ (6) 

In the presence of vortices, ip has singularities at the locations r 4 of each vortex 
(of vorticity Qi = ±1): 

VxV^r)=^2^( 2 '(r-4 (7) 

i 

A natural way to proceed with solving Eq. (1) is to attempt to include the effects of 
vortex excitations perturbatively by expanding in small phase gradients. As ip ap- 
pears only in the exponential of the off-diagonal terms, one is motivated to perform 
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the following naive gauge transformation 



u n —> u n = e lip/2 u n , (8) 



V n -> V n = C tv/2 V„ 



This effectively moves the phase from the off-diagonal components and introduces 
it as a phase gradient in the diagonal components of the BdG equation: 

H'( i ! n )=E n ( i t n ), (9) 



H'= "V W " n i^ a , ..)■ (10) 



-(V + fV^) 2 -^ A 

However, as discussed by Balents et al 10 , such a gauge transformation produces 
branch cuts in the Bogoliubov-de Gennes eigenstates. This can be seen by noting 
that, if we assume that u n and v n are single-valued functions then Eq. (9) must be 
solved under the condition that u n and v n each gain a factor of e 47r upon encircling 
every vortex. The authors of Ref. 10 argue that such branch cuts lead to frustration 
effects that would strongly favor the pairing of vortices. Thus, by only allowing 
doubly quantized vortices, there are no branch cuts in the fields u n and v n to worry 
about. 

Here, we follow Franz and Tesanovic 11 in assuming that the correct approach is 
to consider the condensation of singly-quantized vortices in a d-wave superconduc- 
tor. They avoid the branch-cut difficulty by making a gauge transformation which 
does not directly introduce branch cuts but still keeps track of their physical effects. 
Briefly, their technique 14 involves splitting the vorticies into two groups labelled 
"A" and "B" and then performing a gauge transformation of the form 

u n -» u n e^ A (11) 
v n ^v n e-^ B , (12) 

where <Pa(b) is the phase associated only with the vortices in group A(B). This 
gauge transformation has the advantage of not introducing any branch cuts in the 
quasiparticle wave function while at the same time treating the electrons and holes 
on an equal footing. The Berry gauge field a emerges, upon averaging over all vortex 
configurations, as the difference of phase gradients Vy^ — V<^b. 



3. Path-integral treatment of the Aharonov-Bohm effect 

In Sec. 2 we discussed one difficulty associated with solving the BdG equation in the 
presence of vortex excitations, i.e., that naive gauge transformations of the form of 
Eq. (8) introduce branch cuts in the quasiparticle wavefunctions. To motivate the 
path integral technique which we shall use to handle these branch cuts, in the 
present section we review previously known results 19 > 2 o,2i,22,23 m a re l a teci system 
in which gauge transformations must be made with care: The Aharonov-Bohm 
problem 15 . Our aim is to find an expression for the propagator G(x, y, t) for an 
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electron constrained to a ring; we represent this physical space by the symbol M. A 
thin solenoid penetrates the origin and is described by a vector potential A r = and 
Ag = (f>/2Trr, with <f> being the total flux. The Schrodinger equation for G(x,y,t) 
takes the form 

(-id t - (V - zeA) 2 )G(x, y, t) = 5(x - y)S(t). (13) 

From Feynman's path-integral representation of quantum mechanics, we know that 
G(x, y, t) may be expressed in terms of a sum over all paths of an amplitude for 
each path. In the present situation, the doubly-connected nature of M indicates 
that such paths may be divided into homotopically distinct classes depending on 
how many times the path in a particular homotopy class winds around the origin. 
Thus, following Schulman 19 , we may express G in the form 

G(x,y,i) =]TG„(x„,y,t), (14) 

n 

where the function G„(x„, y, t) contains only Feynman paths that wind around the 
origin n times. The subscript n on x reminds us that x„ is the same as x after 
having wound n times around the origin. Let us consider the physical meaning of 
G„(x n ,y, t). The space M on which Eq. (13) is to be solved has the topology of a 
torus. However, by keeping only Feynman paths which wind n times, it is as if we 
are solving the same equation on a helix-shaped space which winds n times: A helix 
(topologically, a line) is the universal covering space of a torus. We shall denote the 
universal covering space by the symbol M*; the important properties of M* are 
as follows: 1) M* is locally equivalent to M, and 2) M* is simply connected. It is 
further true that the Schrodinger equation [i.e. Eq. (13)] is satisfied for each of the 
terms G n (x n ,y, t) entering Eq. (14) 24 . 

The fact that M* is simply connected means that the quantity 
exp (ie f* n A • dsj is well-defined on it and can be used to simplify the equation 
for G„. Thus, by writing 

G„(x„,y,t) = G„(x„,y,£)exp (ie J A ■ ds^j , (15) 

it can be seen that G„ satisfies 

(-id t - V 2 )G„(x„, y, t) = <5(x„ - y)S(t), (16) 

i.e., it is the Green function for the free-particle (by this we mean A = 0) 
Schrodinger equation on M*, which we denote by G^ (x„, y, t) and has the ex- 
plicit form 20,21,22,23 

G(°)(x„,y,t) = ^ ix2+y2) J JXe^'- e+ ^h w (g) , (17) 

where I\ is the modified Bessel function and we have written the coordinates x„ 
and y in terms of the radial coordinates (x,9' + 27m) and (y,9), respectively. We 
emphasize that this is not the usual two-dimensional free particle propagator, but 
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that part of the propagator due to Feynman paths with winding number n (or, 
equivalently, the free propagator on M*). The usual free-particle propagator may 
be obtained from Eq. (17) by summing over all winding numbers n. 

The next step is to insert G„ (x„,y,i) into Eq. (14) after having directly cal- 
culated the factor exp (ie J* n A • dsj : 

G(x, y, t) = e^ 9 '- e+2 ^^G^ (x„, y, t) (18) 

n 

= 4^*" W 'I>"*'- 8, '|™-.,(1). <»> 

m 

where a = e<fi/2ir measures the flux through the hole. Equation (19) is the single- 
particle propagator for the Aharonov-Bohm problem 20 ' 21 ' 23 j which is obtained from 
Eq. (18) using Eq. (17) along with the Poisson summation formula. 

For our purposes, we are primarily interested in Eq. (18), which exhibits the 
structure proposed by Schulman 19 : The propagator is expressed as a sum of terms, 
each of which is a free propagator associated with a particular winding number 
(i.e. on the covering space) multiplied by a gauge-field dependent factor. Before 
proceeding, however, let us briefly note one feature of Eq. (19), namely that for 
integral a, the summation over m may be shifted by a and Eq. (19) becomes 
the two-dimensional free propagator, i.e., the solution to Eq. (16) in free space. 
This is merely the statement that for integral a there are no interference effects 
associated with the flux and the introduction of the universal covering space was 
unneccesary. For non-integral a, however, the universal covering space provides a 
natural way to identify Feynman paths in different homotopy classes which have 
different phase factors arising from the gauge field. In a similar fashion, we shall 
see that our technique for handling vortex excitations in d-wave superconductors 
by invoking the covering space is useful for the case of hc/2e vortices [for which 
the gauge transformation in Eq. (8) produced branch cuts] but unnecessary for 
doubly quantized vortices [for which the gauge transformation in Eq. (8) produced 
no branch cuts]. In the next section, we shall exploit the universal covering space 
to gain insight into the problem of hc/2e vortices in d-wave superconductors. 



4. Vortices in a d-wave superconductor: Static case 

Having discussed the covering-space approach to the Aharonov-Bohm problem, in 
the present section we apply these ideas to the context of multiple vortices in a 
d-wave superconductor. In this section we shall assume a high-temperature approx- 
imation in which the vortices can be taken to be static. The extension to quantum- 
vortex fluctuations will be straightforward and discussed in Sec. 6. Our starting 
point is the following expression for the BdG propagator averaged over vortex po- 
sitions: 

[G(r,r';r)] v = / V V G(r, r'; r)e' 3 ^, (20) 
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(d T + H)G(r, r'; t) = 5^ (r - r')tf(r), (21) 

where G(r, r';r) is the BdG Green function in the presence of a fixed pattern of 
vortices and S[<p] is the Boltzmann weight associated with a particular, static, pat- 
tern of vortices. The Hamiltonian H is given by Eq. (2); however, henceforth wc 
shall not display the physical gauge field A. It may be reinstated at the end of the 
calculation. 

Each vortex configuration entering into Eq. (20) corresponds to a particular 
set of locations for the in Eq. (7). Thus, as in the Aharonov-Bohm problem, 
the manifold M on which we aim to solve Eq. (21) is multiply connected, i.e., 
it has "holes" at each r^. To ascertain the physical importance of this multiple- 
connectedness, we exchange the original multiply-connected space M (on which 
Eq. (21) is defined) for its simply-connected covering space M* on which a gauge 
transformation of the form of Eq. (8) may be performed. As we saw in Sec. 2, 
for the Aharonov-Bohm effect the relevant universal covering space was a line. In 
the present case, the fact that we shall always consider the case of multiple vortices 
means that the associated universal covering space is difficult to envision physically. 
Fortunately, as in the case of the Aharonov-Bohm effect, the covering space will 
merely serve as a way to keep track of the phase picked up by various Feynman 
paths. 

Following our procedure in the preceding section, we express G(r, r'; r) as a sum 
of propagators, the Feynman paths of which are grouped into topologically distinct 
sectors: 

G(r,r / ,r) = ^G n (r n ,r / ;r), (22) 

where the vector index n = (ni,n 2 ,---) keeps track of the winding number rn 
associated with the ith vortex. Technically speaking, this winding number does 
not fully label the elements of the non-abelian homotopy group associated with the 
multiply-connected space in question, but this does not matter for the calculation we 
are attempting since all that matters is that there exists a label for these elements. 
On M* we have at our disposal the quantity 

Xn(r„)= / ° Vp(r) • dv, (23) 
Jo 

which we emphasize is not uniquely defined on the space M. For simplicity, we have 
set r' = 0. Physically, %n represents a line integral along a particular path through 
the vortices. We write G n (r n ,0,r) in the form 



G n (r n ,0,r) = exp 



^Xn(r n )<r 3 



G„(r n ,0,r), (24) 



where 03 is the usual Pauli matrix. As in Sec. 3, this amounts to making different 
gauge transformations for Feynman paths having different winding numbers. The 
G n satisfy 

(d T + H')G n (r n ,0:T)=5^(r n )5(T), (25) 
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Fig. 1. Sketch of various homotopically incquivalcnt Fcynman paths (black curves) connecting 
the points y and x in the presence of vortices (represented by large crosses). As an example of 
how the winding number n is defined, we note that for the leftmost vortex n can be taken to be 
for path A and 1 for path B. 



Having made a gauge transformation on the covering space, we proceed by com- 
puting the path-dependent gauge-transformation factor exp [|xn(r n )o'3] . The path- 
dependence of Xn(r n ) is exhibited schematically in Fig. 1. Consider for example the 
loop made by paths B and C: since there is one vortex inside the loop, we have 



/ ■ 

Jb-c 



Vip{r) ■ dr = ±27r, 



(26) 



depending on the vorticity of the enclosed vortex. Thus, Xn( r n) differs by 2w for 
paths B and C. More generally, we may define 

Xn(r n ) = <^(r)-^(0) + 2^^ ft n J , (27) 



5Xn(r„)CT 3 _ e ^(¥>(r)-y(0))*3 p i7r£i qiTH 



e 2 ' 



(28) 

where it is important to note that in Eq. (28) we have used the fact that 
cxp(i7m(73) = cxp(i7rn) for integer n. We thus have the following expression for 
the averaged BdG Green function: 



[G(r,r';r)] v 



/ 



Vipe 



S[v] e 5V( r )*3 



^G„(r n ,r';r)e-S 



e 2 



|¥>(r')d- 3 



• (29) 



Since we are considering only singly quantized vortices, so that qi = ±1, the branch- 
cut factor expi7r^ i qirii is given by ±1 for different trajectories and is the manifes- 
tation, within the present approach, of the quasiparticle branch cuts encountered in 
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Sec. 2. In Fig. 1, since they are separated by two vortices, the branch-cut factor is 
the same for paths A and B. In general, if we considered the condensation of pairs of 
vortices (g, = ±2), then these factors could be all taken to be unity and there would 
be no need to invoke the covering space construction. However, although there is 
a Z 2 character to the physics of fluctuating branch cuts as discussed here, it has a 
origin which is distinct from that of the Z2 gauge theory of Senthil and Fisher 25 , 
which, like the theory of Ref. 10, only incorporates doubly quantized vortices. 

Before proceeding to the next stage of the calculation, we pause to note that 
the argument of the functional integral in Eq. (29) is the analogue, within the 
present context, of Eq. (18) in the calculation of the Aharonov-Bohm propagator 
in Sec. 3. For that case it is possible to calculate the propagator on the covering 
space [cf. Eq. (17)] explicitly. In the present case, however, G„(r n ,r';r) is most 
likely not analytically solvable and, in the next section, we proceed by making some 
simplifying approximations. 



5. Effect of branch cuts on quasiparticle dynamics 

In the present section, we shall attempt to evaluate the functional integral over 
vortex and spin- wave like phase fluctuations in Eq. (29). As discussed in Sec. 4, the 
effect of the branch cuts in the quasiparticle wavefunctions has been traced back to 
the "branch-cut" factor expi-Tr^ ^rij in Eq. (28) which differs for homotopically 
inequivalent trajectories. If this factor were absent (e.g. if we had only allowed 
the proliferation of doubly-quantized vortices so that, in effect, J2i 1i n i would be 
constrained to be an even integer), then there would be no need to divide the 
equation for the propagator into topologically inequivalent sectors, and the gauge 
transformation given in Eq. (8) could have been made on M. 

Let us turn to the evaluation of the functional integral over phase fluctuations. 
Formally, ip appears in two distinct places in Eq. (29): Firstly, it appears in the 
definition of H' in Eq. (10). As our aim is to focus on the effect of branch cuts, here 
we neglect V<^ in H' . This simplifying approximation is not technically necessary at 
this point. It is motivated by the fact that, were we to keep the local phase gradient 
as another ("Doppler") gauge field v then, as discussed in Refs. 11, 12, v would 
end up being irrelevant in the rcnormalization group sense. Secondly, <p appears in 
the exponential factors exp (±i(p(r)as/2) in Eq. (29). Let us examine the matrix 
structure of the argument of the functional integral by writing it as a 2 x 2 matrix: 

i , - i 1 ( G„ lie 2M r )-v( r ')) q ^eKvW+vCr')) 



E, , ( Gr, netwW^t 1 ' )) Gr, i9e^ (v(r)+¥,(r )•> \ 



By writing out the full matrix form for the product of these three factors appearing 
in Eq. (29) we sec that whereas the diagonal components contain a phase difference, 
the off-diagonal components contain a sum of phases and thus must vanish upon 
performing the average over the <p. This merely implies that the off-diagonal com- 
ponents of the single-particle Green function must vanish in the pseudogap regime, 
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i.e., that superconductivity has been destroyed. Being phase differences, the phase 
factors appearing in the diagonal components do not vanish upon averaging over ip. 
We shall take them to be unity, invoking the same argument which we used for the 
phase gradients in H'. 

By making the preceding approximations, we are assuming that the branch cuts 
factors lead to quantum interference effects that have a dominant effect on the 
low-energy quasiparticles. We thus arrive at 



[G(r,r';r)]^ f 2tye- s M£\ 



,"r£ 4 9.". 



x (r n ,r 



8 T + 



-V 2 -^ A p x p y 

AoPxPy V 2 + fl 



|r',0) 



(31) 



diag. 



where we have explicitly written the expression for G„(r n , r'; r) within this approx- 
imation. The subscript "diag." indicates that our expression for [G(r, r'; t)] v only 
includes the diagonal components of the expression on the right side of Eq. (31), the 
off-diagonal components being zero as noted above. Henceforth, we shall suppress 
this subscript, although all subsequent equations share this property that only the 
diagonal components are of physical interest. 

The next step is to represent the branch-cut factor in terms of the flux of a 
field a, allowing us to evaluate the functional integral over vortex configurations in 
Eq. (31). To do this, we define a such that its local "magnetic field" has singularities 
at the locations of all the vortices (but curl- free on M*): 



V x a(r) = 7r^g i( 5 (2) (r-r;). 



(32) 



To account for the factor e I7r ^ < qiUi ,we write qirii as the flux through a carefully 
chosen loop: 



qirii = -— ( a • dr + — [ a • dr, 
V n Jn fir 



(33) 



where the subscript n indicates that the first integral winds along the path falling 
in the homotopy class labelled by n (i.e. from y to x„ in M*), and the subscript T 
indicates that the integral is to be taken along some arbitrary fixed reference path 
xr (a): 

/ a • dr = / dsa fl (x r {s)) : 
Jr Jo 

where x r (0) = r' and x r (l) = r. Thus, we have 



ds 



(34) 



[G(tJ;t)] v 
x(r„,r| 



J V^e~ s ^ J 



Vac 



if r a-d 



-V 2 -^ A P X Py 

AqPx Py V 2 + /I 



i J n 



k',o) 



■dr 
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x^VxaW-i^^'tr-r^ (35) 

i 

where we must emphasize that the field a has entered in a gauge-invariant way, in 
the sense that [G(r, r'; t)] v is manifestly dependent only on V x a. Although the 
quantity / a - dr in Eq. (33) is trajectory-dependent, the integrals along the path T 
are by definition the same for each term in Eq. (35). Before proceeding, we remark 
that this propagator is reminiscent of the gauge-invariant Green function discussed 
in Ref. (16). 

Next, we make another gauge transformation of the form of Eq. (24), absorbing 
the factors exp —i J n a • dr into the Green function associated with that particular 
trajectory (or, rather, class of trajectories associated with a particular homotopy 
class in the presence of the vortex excitations) . After making such a gauge transfor- 
mation, the division of propagators into topologically distinct sectors is redundant 
(the field a keeps track of the branch cut factors), leaving us with the following 
expression for the Green function: 

[G(tJ;t)] v * J !><pe- s Mj Pae^r »^ (r , T | ^ + ^ | r ', ) 

^(VxaW-^^^r-ri)), (36) 



H[*]=l^*J-» (37) 



(p + a) 2 - fjL D 

I) -(p + a) 2 +/i 

D = ^-[(Px + a x )(Py + Cly) + (p y + Oy ) (j) x + (l X )\. (38) 

The most straighforward way to verify the equality of Eq. (35) and Eq. (36) is 
to examine the gauge-transformation properties of the arguments of the respective 
functional integrals. The next step is to evaluate the functional integral over vortex 
and spin- wave like excitations, leaving an effective functional integral over the field 
a. Formally, we may follow Franz and collaborators 11 , writing 



J V V c- s ^5{V x &{v)-2n^q l 5 {2 \v-Y l )). 



(39) 



Using this definition of the vortex action, we arrive at our final expression for the 
single-particle BdG propagator in the presence of static vortex excitations: 

[G{v,v'-r% ~ J Pac- S *e* W (r , T |_^__| r ', >. (40) 

Equation (40) is essentially our main result, i.e., we have obtained a gauge-theory 
model for vortex excitations in a d-wave superconductor via a path-integral tech- 
nique. However, we recall that until now we have considered only static vortex 
excitations, so that the functional integral over vortex excitations in, e.g., Eq. (20), 
is truly a thermal average. As we shall see, however, the generalization to arbitrary 
dynamic vortex excitations is straighforward within the path-integral technique. In 
the next section, we discuss this issue in detail. 
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6. Vortices in a d-wave superconductor: Dynamic case 

In arriving at Eq. (40), we have made use of the Feynman path integral description 
of the BdG propagator in the presence of a static pattern of vortices. In the present 
section, we generalize this procedure to the case of arbitrary fluctuating vortices. 
The quantity we are interested in calculating is, formally the same as was discussed 
in Sees. 4 and 5 and is given by 

[G(r,r';r)] v ee J V V c' s ^ (r, r|^-^ |r', 0), (41) 
-V 2 - \i e ^/ 2 Ae^/ 2 ' 

where the only distinction from the preceding discussion is that now tp is time- 
dependent and S[(f\ is an action, as opposed to a Boltzmann. 

In the picture discussed in the preceding sections, the various Feynman paths 
were parametrized by the temporal variable and fell into homotopically incquivalent 
classes defined by the locations of the (static) vortices. Now imagine that the vor- 
tices are time-dependent. One may of course still construct a conventional Feynman 
path integral even though H is time-dependent. However, since the vortex positions 
evolve in time, Feynman paths parametrized by time do not fall into homotopically 
inequivalent classes in the same way. Thus, it is convenient to use a path-integral 
method which treats the time variable r on an equal footing with the spatial coor- 
dinates. This may be realized by using a so-called "fifth-parameter" path integral 
scheme 19 . Within such a scheme, one introduces a ficticious temporal coordinate 
A and constructs a path integral representation of the corresponding Green func- 
tion in which the temporal variable in the Feyman paths is A, not r. For example, 
consider the following Schrodingcr equation: 

(id x + C) T(x, y- A) = -iS(x - y)S(X), (43) 
C = 8 T + H, (44) 

where x = (r, r) and y = (r',0). It is clear that the propagator F(x,y;\) may be 
expressed as a path integral in which the paths are labeled by A. Furthermore, it 
may be simply related to the argument of Eq. (41) via 

1 r°° 

( r > T \jT—Fj\ r '°) = / d\r(x,y;\)e- e \ (45) 



o 



1 d T + H 1 
where e = 0+. 

The utility of the preceding discussion is that now we have a path-integral rep- 
resentation (i.e. that of Eq. (43) for the argument of Eq. (41) in which the variable 
r is treated like any other coordinate, so that instead of thinking of vortices in 
a two-dimensional XY model we envision vortex loops in a three-dimensional XY 
model. The vortices are "static" relative to the coordinate A, allowing us to split the 
Feynman paths into homotopically distinct sectors depending on how they wind rel- 
ative to them. This amounts to constructing an equation for T(x, y; A) that is of the 
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form of Eq. (22). Following the same steps as in the static case (i.e, making gauge 
transformations on the covering space, etc.), generalized from 2 to 3 dimensions, wc 
find that [G(r, r'; t)] v is given by 

[G(r, r';r)]^J Va e~W Jr -*(,, r\ ^ + ^ + ^ \r>, 0) (46) 

S a = -K{d v a» - (47) 

where now the vector field a has three components given by a — (a T ,a x ,a y ). The 
boldface quantity a still refers to the spatial vector a = (a x ,a y ), and H[a] is still 
given by Eq. (36) (although we again emphasize that now a fluctuates dynamically). 
The line integral with subscript label T is similarly generalized from Eq. (34) to now 
indicate a line integral from (r',0) to (r, r) The gauge-invariant action S a (with K 
being an appropriate coupling) is formally defined by an equation exactly analagous 
to Eq. (39); in Eq. (47) we have used the results of Franz et al. 11 who computed it 
using (see also Ref. 12) a model of free fluctuating vortex loops. It is important to 
note however that this form for S a could be expected on general symmetry grounds, 
it being the simplest quadratic action which respects gauge invariance. 

To make the connection between [G(r, r'; t)] v and the propagator for the QED 3 
gauge theory (i.e. to make the gauge symmetry manifest) we represent [G(r, r'; t)] v 
in terms of a fermionic path integral 

[G(t,t';t)] v ~ j V^V^Va^ s ^^ s ^{Y,T)^{r l ,Q)e t ^^ d ^ (48) 

Sip= J d 3 x^ [8 T + ia T + H[a}} tp, (49) 

where tp^ and tp represent Nambu spinor fields. The action S,p + S a is invariant 
under the gauge transformation ip — > c~ lx ip, — > a M + d^x- I n addition, we see 
that the factor e* ^r a dr in the argument of Eq. (48) ensures that 

V'(r,T)V t (r',0)e I .''r a - dr , (50) 

and therefore [G(r, r'; r)] v , is also gauge- invariant. One may also expand the 
fermions near the nodes of the d-wave order parameter; by combining the Nambu 
fields ip m to four component Dirac spinors 11 one may complete the connection to 
QED 3 . However, for the present purposes it is sufficient to have the above gauge 
theory be our final expression. 

7. Discussion 

In this Paper, we have explored the effect of vortices on the the quasiparticle exci- 
tations of d-wave superconductors in an attempt to understand the pseudogap phe- 
nomena. Our final expression for the single-particle Green function (i.e. Eq. (46)) 
makes direct contact with the QED 3 11 scenario of the pseudogap phenomena. As 
discussed in Sec. 1, the original formulation of the QED 3 of the pseudogap phe- 
nomena 11 relied on splitting the vortices into two groups (A and B) and making a 
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gauge transformation of the electrons with respect to group A and the holes with 
respect to B. The aim of such an approach is to avoid directly introducing branch 
cuts in the quasiparticle wavcfunctions. 

The approach presented here deals with the quasiparticle branch cuts in a more 
direct fashion, by realizing that the problem of solving the BdG equation in the 
presence of vortices may be formulated as a problem of making gauge transforma- 
tions on multiply connected manifolds. By going to the universal covering space of 
the associated manifold, the branch cuts emerge in a controllable fashion as factors 
of ±1 multiplying various Feynman paths. The field a arises in an attempt to mimic 
the effect of these branch-cuts; in particular the line integral factor in Eq. (48) comes 
from expressing the branch cut factors ±1 associated with various Feynman paths in 
terms of the flux of a. Although we have clarified the origin of the Berry gauge field 
in the BdG action, we have not clarified issues associated with the calculation of 
the single particle Green function [i.e. Eq. (48)] As has been noted in several recent 
papers 18 > 26 > 27 7 there are many technical difficulties associated with the calculation 
of this quantity due to the presence of the line integral over the gauge field. As 
noted in a recent preprint 28 , however, for certain gauge- invariant quantities these 
line integrals do not appear. 
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